void LowPass (Rad, ImSize, Image)
int   Rad, ImSize;
float ** Image;
{
   int  r,c, Bot,Up;
   long rr,cc;

   Bot = ImSize/2 -1;
   Up  = ImSize/2 +1;
   for (r=-Bot; r<Up; r++)
      for (c=-Bot; c<Up; c++)
        if (Rad < (int) sqrt ((double) r*r+c*c))
           Image [r+Bot] [c+Bot] = (float)0;
}

void HighPass (Rad, ImSize, Image)
int   Rad, ImSize;
float ** Image;
{
   int  r,c, Bot,Up;
   long rr,cc;

   Bot = ImSize/2 -1;
   Up  = ImSize/2 +1;
   for (r=-Bot; r<Up; r++)
      for (c=-Bot; c<Up; c++)
         if (Rad > (int) sqrt ((double) r*r+c*c))
            Image [r+Bot] [c+Bot] = (float)0;
}
